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The scatter in the radial profiles of X-ray luminous galaxy clusters 
as diagnostic of the thermodynamical state of the ICM 
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ABSTRACT 

We study the azimuthal scatter in the radial profiles of X-ray luminous galaxy clusters, with 
two sets of high-resolution cosmological re-simulations obtained with the codes ENZO and 
GADGET2. The average gas profiles are computed for different angular sectors of the cluster 
projected volume, and compared with the mean cluster profiles at each radius from the center. 
We report that in general the level of azimuthal scatter is found to be ~ 10 per cent for gas 
density, temperature and entropy inside R200, an d ~ 25 per cent for X-ray luminosity for 
the same volume. These values generally doubles going to 2i? 2 oo from the cluster center, and 
are generally found to be higher (by ^ 20 — 40 per cent) in the case of perturbed systems. A 
comparison with results from recent Suzaku observations is discussed, showing the possibility 
to simply interpret the large azimuthal scatter of observables in light of our simulated results. 
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Only a fraction of the volume occupied by the Intra- 
Cluster Medium is currently mapped with X-ray obser- 
vations. The typical measures of the surface brightness, 
temperature and metal abundance extend out to a frac- 
tion of the virial radius, R v i r , and only in a few cases 
meaningful estimates at -R500 B or beyond are obtained. 
Because of this, more than two-thirds of the typical 
cluster volume, where the primordial gas is accreted 
and the dark matter (DM) halo is forming, are still un- 
known for what concerns the gas/DM mass distribution 
and the related thermodynamical properties. This poses 
a limitation in our ability to characterize the physical 
processes presiding over the formation and evolution 
of clusters, and thus to use them as high-precision cos- 
mological tools (see Ettori & Molendi 2010 and refer- 
ences therein). 
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1 ^500 (R200) i s the radius enclosing a mean matter over-density of 500 
(200) times the critical density of the Universe. Compared to the cluster 



virial radius these radii correspond to -R500 
0.7-Rvir. 



0.5-R v i r and to -R200 



In general the surface brightness distribution, re- 
quiring lower counts statistic, is easier to be charac- 
terized and has been successfully studied with Rosat 
PSPC (e.g. the outskirts of the Perseus cluster in Et- 
tori et al. 1998; more systematic studies have been car- 
ried out from Vikhlinin et al. 1999 and Neumann 2005) 
and with Chandra (e.g. Ettori & Balestra 2009). On 
the other hand the measure of the gas temperature is 
more strongly limited from the background emission. 
Current exposures with Chandra (e.g. Vikhlinin et al. 
2006) and XMM (e.g. Leccardi & Molendi 2008) en- 
abled the reconstruction of the azimuthally-averaged 
temperature profile for few tens of nearby galaxy clus- 
ters out to ~ i?5oo- A step forward has been recently 
obtained from the observations of a handful sample of 
nearby X-ray luminous objects with the Japanese satel- 
lite Suzaku that, despite the relatively poor PSF and 
small FOV of its X-ray imaging spectrometer (XIS), 
benefits from the modest background associated to its 
low Earth orbit (see results on PKS0745-191, George 
et al. 2009; A2204, Reiprich et al. 2009; A1795, Bautz 
et al. 2009; A1413, Hoshino et al. 2010, A1689, Kawa- 
harada et al.2010). However as consequence of the 
small FOV (and of the large solid angle covered from 
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Figure 1. Projected X-ray luminosity from each pixel at [0.5 — 2]keV for 
cluster E14 of the ENZO sample at z = 0. The left panel is for the original 
data, the right panel is for the data after the 99 per cent filtering. The side of 
the slices is ~ 3 X 5Mpc/h, the pixel size is 25kpc/h. 



the bright nearby clusters) these observations provide 
often discrepant constraints over a few arbitrarily cho- 
sen directions. 

A number of numerical and semi-analitical works 
in the last few years aimed at providing reliable mod- 
els of the thermal properties of galaxy clusters at their 
outskirts, from the theoretical point of view (e.g.Loken 
et al.2002; Roncarelli et al.2006; Burns, Skillman 
& O'Shea 2010; Lapi, Fusco-Femiano & Cavaliere 
2010; Cavaliere, Lapi & Fusco-Femiano 2010; Vazza 
et al.2010). Quite recently, a detailed analysis of the 
cluster profiles up to « 2 • R200 was reported by Burns 
et al.(2010), using a sample of clusters simulated with 
the ENZO code (see Sec©. 

In this article, we aim at characterizing the intrinsic 
azimuthal scatter (due to the 3-D structure of clusters) 
that should affect at some level the measurements of 
the thermodynamical properties in the cluster outskirts, 
and at providing a guidance in interpreting the existing 
and incoming observations. To this end, we extended a 
method proposed by Roncarelli et al.(2006) to compute 
the cluster profiles by avoiding the bias caused by small 
collapsed gas clumps, and we measured the level of 
azimuthal scatter around the average clusters profiles 
in a total sample of 29 massive galaxy clusters. 

The paper is organized as follows: in Secf2]we de- 
scribe the set of simulations; in Secj3]we present the 
average properties of the azimuthal scatter of the gas 
density, temperature, entropy and X-ray luminosity, 
with a comparison with the observational constraints 
available. The discussion of our results is reported in 
Secf4] In the Appendix, additional test on the physical 
and numerical convergence of the main results of the 
paper are discussed. 
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Figure 2. Filtered profiles of gas density (top), gas entropy (center) and gas 
temperature (bottom) for all GADGET clusters (blue) and for all ENZO 
clusters (red). The profiles of all quantities are normalized at i?200- The 
dotted/dashed lines refer to each cluster in the sample, the thick solid lines 
show the average profiles for each of the two samples. 
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Table 1. Numerical and cosmological parameters for the simulations: code, 
mass resolution for DM, gravitational softening, maximum gas spatial res- 
olution <5 max , cosmological DM density, cosmological baryon density, nor- 
malization of the primordial density power spectrum. GADGET2 runs 
adopted h = H/(100km/(s ■ Mpc)) = 0.7 while ENZO runs adopted 
h = 0.72. 



code 



[M Q /h] [kpc/h] [kpc/h] 



ENZO 6.7- 
GADGET2 1.13 



10 s 
• 10 9 



50 
5 



25 
- 10 



0.226 
0.261 



0.044 
0.039 



2 NUMERICAL METHODS 



2.1 Cluster simulations 



We analyze two independent sets of high-resolution 
(<5max ~ 10 — 20kpc/h) cosmological simulations 
of galaxy clusters, extracted from large cosmologi- 
cal boxes (Lbox ~ 500 Mpc) and re-simulated with 
two of the most tested numerical codes for high per- 
formance computing: the Adapative Mesh Refinement 
code ENZO and the Tree Smoothed Particles Hydrody- 
namics code GADGET2. The main cosmological and 
numerical parameters adopted are listed in Table 1. 

A data-set of 20 non-radiative re-simulations of 
galaxy clusters with total masses 6 • 10 14 ^ M/M & ^ 

2 • 10 15 was produced with the 1.5 public version of 
ENZO 1.50(Norman et al.2007 and references therein), 
with a few physical implementation (Vazza et al.2010, 
hereafter ValO). For the specific details about the 
simulational setup and the post-processing procedure 
adopted we address the reader to ValO; a public archive 
of part of this data-sample is currently available at 
http://data. cinecaJi\ 

Massive halos in ENZO runs were identified fol- 
lowing a spherical over-density method, applied to the 
simulated outputs at z = (e.g. Gheller, Pantano & 
Moscardini 1998). All clusters with total masses larger 
than ~ 6 • 1O 14 M were extracted and sampled at the 
highest available resolution in the simulation 25kpc/h 
out to the distance of ~ 3 — 4i?2oo from the cluster 
centers. 

A second sample of 9 non-radiative re-simulations 
of galaxy clusters with masses 2 • 10 14 ^ M/Mq ^ 

3 • 10 15 was obtained with GADGET2 (Springel 2005). 
The details for the simulational setup of these runs 
can be found in Dolag et al.(2005) and Roncarelli et 
al.(2006), hereafter Ro06. 

Also in this case, the halos at z = were identi- 
fied basing on a spherical over-density method; the ex- 
traction procedure was tailored in order to have a sam- 
ple of obviously non-merging systems across the mass 
range, discarding double systems or systems with on- 
going strong interactions. 



2 ENZO is developed by the Laboratory for Computational Astrophysics 
at the University of California in San Diego (http://lca.ucsd.edu I 
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Figure 3. 3-D profiles of gas density for one perturbed ENZO cluster (top) 
and one relaxed ENZO cluster (bottom). Different colors represent profiles 
for 8 different sectors, the solid lines are for the un-filtered cells while the 
dotted lines are for the 99 per cent filtering. The red line in overlay shows 
the best-fit profile for the filtered data of all GADGET2 clusters (from in 
Roncarelli et al.2006). 



2.2 Data reduction. 

In the present work, we study the distribution of the 
ICM density, temperature, entropy and surface bright- 
ness in our two sets of simulated galaxy clusters by 
analyzing their radial profiles in the range between 
0.01 ^ r/R2oo ^ 2. For both set of simulations, the 
radial profiles were computed within shells of equal 
spacing (rj 25kpc/h) centered on the center of total 
(gas+DM) mass of each system. The gas matter distri- 
bution in the cluster volume is usually characterized 
by the presence of very dense sub-clumps, with the 
typical dimension of galaxies or small galaxy-groups 
(M < 1O 13 M ). In real X-ray observations, these 
dense and bright clumps can efficiently be masked out 
when present, and since in this work we aim at com- 
paring with observational results, we applied a recipe 
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Figure 4. Scatter profiles of projected Lx at [0.5-2]keV for two represen- 
tative clusters of sample (ENZO=top, GADGET2=bottom). The solid lines 
are for the un-filtered cells while the dotted lines are for the 99 per cent 
filtering. Every color refers to a different number N of angular sectors in 
the projected volume. 

to mimic the excision of such clumps from the simula- 
tions too. 

We followed the procedure introduced in Ro06 and 
for every 3-D shell in the cluster radial profile we 
sorted the gas particles/gas density in decreasing vol- 
ume order (thus in increasing density order). We next 
computed the profiles of gas density, of spectroscopic - 
like temperature (Tsl, see Rasia et al.2005), of entropy 
and of X-ray luminosity in the [0.5-2]keV band con- 
sidering only the particles/cells summing up to a fix 
percentage of the total volume of the shell. For con- 
sistency with Ro06, this percentage was set to the 99 
per cent of the total volume of each shell. This essen- 
tially make sure that, at each radius, the first densest 
percentile of the gas distribution within each spheri- 
cal shell is removed from the data-set, and replaced 
with the average value at the same radius. Tests showed 
that the excision of this first percentile from the origi- 
nal data-set efficiently removes most of collapsed gas 
clumps at each radius. For further details we address 
the reader to the discussion in Ro06. A similar filtering 



procedure was applied also directly to the 2-D maps of 
projected X-ray luminosities, by substituting the 1 per 
cent most luminous pixels in all maps with the aver- 
age value at the corresponding radius. Figure Q] shows 
the application of this method for one projected X- 
ray map of one cluster from the ENZO sample; sim- 
ilar maps for the GADGET clusters can be found in 
Ro06. In Figure [2] we report the gas density, TgL and 

gas entropy (denned as S = Tg^/n^ 3 , where n e is 
the electron number density) for all clusters in the two 
samples, with the 99 per cent filtering (all profiles are 
normalized at i?200 to highlight self-similarity). Since 
the ENZO sample contains more clusters, larger sta- 
tistical fluctuations are found in this data-set; neverthe- 
less the average profiles obtained for the two numerical 
methods are in agreement within a 10 — 20 percent for 
r > 0.2 ■ i?200- On the other hand, as expected siz- 
able differences are found for the innermost cluster re- 
gions, r < 0.2 • i?200> in li ne with the findings reported 
and investigated by many different works in the liter- 
ature (e.g. Frenk et al.1999; Voit et al.2005; Wadsley 
et al.2008; Mitchell et al.2009; Springel 2010; Vazza 
2010). Discussing the above differences is far beyond 
the purpose of this work; we just notice that for our 
main goal here the most interesting regions under study 
are > 0.5 • i?200> where a general good consistency be- 
tween the two codes is ensured. The effect of additional 
physical processes not modeled in these runs (e.g. cool- 
ing, star formation, magnetic fields etc.) is not expected 
to play an important role. For completeness, the issue 
of non-gravitational physics in the simulations is in- 
vestigated for a few representative clusters in the Ap- 
pendix (see also Ro06 and discussion therein). 

The simulated X-ray emission from our set of clus- 
ters was produced by adopting a MEKAL emission 
model (e.g. Liedahl et al. 1995) for each gas parti- 
cle/cell, according to its temperature and by assuming 
a constant gas metallicity of Z/Zq = 0.2. This value 
is likely lower than the actual metallicity in the inner- 
most clusters region (e.g. Leccardi & Molendi 2008), 
however it can be considered a good approximation 
to model the spectral energy distribution radiated by 
the ICM at the cluster outskirts, once that the effect of 
metal cooling is considered. 

In order to highlight the systematic effects due to 
the dynamical state of the clusters in our sample, we 
estimated the power ratios P3 / PO, from the multipole 
decomposition of the X-ray surface brightness images, 
versus the centroid shift, w, evaluated within -R500 (see 
Cassano et al.2010). We classified as "disturbed" sys- 
tems those for which the values of P3/P0 > 10~ 7 
and w > 0.02 were found in at least 2 of the 3 pro- 
jected maps along the coordinate axes, or as "relaxed" 
otherwise 0. This classification method splits our two 
samples in the following way: in GADGET2 we ob- 
tain 5 relaxed systems and 4 perturbed systems, while 
in ENZO we obtain 1 1 relaxed systems and 9 perturbed 

3 Since non-radiative runs tend to produce too much concentrated cluster 
cores than real clusters, we expunged the r < O.IR200 region from our 
analysis. 
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systems. Compared to the dynamical classifications in- 
troduced in ValO for the ENZO sample (which were 
based on the matter accretion history inside the virial 
volume, and not on X-ray properties as in this work) 
we report that the "perturbed" class here almost per- 
fectly overlap with the class of "post-merger" systems 
there, while the "relaxed" class contains the "relaxing" 
and "merging" classes of ValO. 

The key feature discussed in this paper is the az- 
imuthal scatter in the cluster profiles along different 
cluster sectors, as a function of the radial distance from 
the center of clusters. We mimicked an observational- 
like procedure and fixed a line of sight towards the 
cluster center, defining a cylinder of radius 2 • -R200 and 
height 4 • -R200 along the line of sight . The volume of 
each cylinder was divided by cutting N "sectors" along 
the line of sight with angular extension A a = 2-k/N . 
This was made in order to approximately reproduce the 
usual decomposition of the FOV in Suzaku observa- 
tions of the peripheral regions of galaxy clusters (e.g. 
Aa ?» 7r/2 in George et al.2009). As an example, we 
report in Fig|3] the behavior of the gas density pro- 
files for a relaxed and a perturbed ENZO cluster, along 
8 different sectors of the projected volume, with and 
without the "99 per cent filtering". For comparison we 
also overplot in within the same Figure the mean gas 
density profile derived in R0O6 after the "99 per cent 
filtering" method (dashed red line). In both cases, the 
profile in at least 2 sectors outside of ~ R v i r shows 
quite large variations, reaching a factor ~ 10 in the 
case of one particular sector in the perturbed cluster. 
These large fluctuations are due to clumps and dense 
filaments crossing the virial radius, and when the "99 
per cent filtering" is applied the agreement between the 
various sectors is generally better than a ~ 20 — 30 
per cent at all radii. A qualitatively similar behavior 
is found also for the other physical quantities, and in 
GADGET2 runs. 

The azimuthal scatter along the cluster radius, 
S c {r), is quantified as: 



Sc(r) 



\ 



N ' 2- 



[Y(r)f 



(1) 



where yi(r) is the radial profile of a given quantity, 
taken inside a given /-sector, and Y(r) is the average 
profile taken from all the cluster volume. In princi- 
ple, functional fitting procedures and ensemble aver- 
ages can be use to derive Y(r) for each cluster, limit- 
ing the bias from substructures in the estimate of the 
average spherical profiles (e.g. Kawahara et al.2008). 
However, we want to make our treatment here the most 
straightforward as possible, in order to be readily use- 
ful to single-object observations, and therefore in this 
work Y{r) is computed directly from the simulated 
cluster cells/particles (or pixels in the case of Lx) af- 
ter the "99 per cent filtering", and with the addition of 

4 In the Appendix we also show a test using the larger line of sight of 
6 • -R200 , showing that our results on the azimuthal scatter of all physical 
quantities are remain unaffected. 



a median smoothing filter along the radius. Our tests 
confirm however that this issue is not very relevant for 
the problem discussed here, and that a simple estimate 
of Y(r) from the filtered data does not affect our statis- 
tical results on S c (r) in any practical way (see SecfJ). 

The choice of the number N of sectors is arbi- 
trary, and may be tuned in order to mimic what done in 
current Suzaku observations. In the following, we will 
show examples of the average behavior of S c (r) using 
N = 2, = 4 and = 8. However, we find that a good 
convergence of S c (r) along the radius, for all analyzed 
quantities, is achieved adopting N ^ 8 for both the dy- 
namical classes considered. As an example, we show 
in Figj4]the azimuthal scatter profiles for the projected 
X-ray projected luminosity at [0.5-2]keV for a cluster 
of each sample. In the "filtered" case, at all radii S c (r) 
do not change by more than a ~ 5 per cent if the num- 
ber of sectors is increased from N = 8 to N = 16. 
In the un-filtered case, however, the presence of the 
point-like bright contribution by collapsed gas clumps 
leave strong fluctuations in the azimuthal scatter even 
at the small angular scale of Aa = tt/8 sampled tak- 
ing N = 16. This may serve as a guidance for future 
X-ray observations of the outermost regions of galaxy 
clusters: an accurate reconstruction of the thermal gas 
distribution of galaxy clusters is reached with good 
accuracy for sectors of angular size Aa ~ 7r/4 (or 
smaller), provided that an efficient masking of point- 
like contribution from dense gas clumps is performed. 



3 RESULTS 

Our total sample of 29 clusters, together with the filter- 
ing procedure outlined in the previous Section, allows 
us to assess the intrinsic level of azimuthal scatter ex- 
pected for the clusters profiles, as an effect of the hier- 
archical clustering and cosmic evolution process driven 
by gravitational processes alone. 

In Figure [5] we report our results for the azimuthal 
scatter profiles of the galaxy clusters simulated in the 
2 codes, by adopting the preferred value of N = 8 
sectors and the "99 per cent filtering" (see Sec l2.2l ). 
Each sample was also divided into perturbed and re- 
laxed systems (see different line-styles in Figj5] Each 
couple of lines with the same style shows the con- 
tours of S c [r) ± o-gr r \ (the standard deviation is com- 
puted across for each dynamical sub-sample of the 
cluster sample). The filtered profiles show an agree- 
ment between the two codes better than a ~ 5 per 
cent at any radius when relaxed systems are consid- 
ered, and slightly larger differences for the case of 
perturbed systems. In general, we report that for gas 
density, gas temperature and gas entropy the intrinsic 
mean azimuthal scatter computed for ./V = 8 sectors is 
S c (r) ~ 10 ± 5 per cent for r < -R200 f° r both codes, 
and S c (r) ~ 25 ± 10 for the X-ray luminosity at [0.5- 
2]keV within the same radius. Outside of R200, both 
codes present a gradual increase of the azimuthal scat- 
ter in all quantities, up to a factor ~ 2 at 2 ■ i?200- In 
general, pertubed objects present larger azimuthal scat- 
ter at all radii, ~ 20 — 40 percent more compared to the 
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Figure 5. Average scatter profiles for all clusters of the paper (ENZO=red; GADGET=blue) for the data with "99 percent filtering" and for TV = 8 sectors. 
For each cluster sample we show the 2 lines enclosing the ±lcr around the mean azimuthal scatter. The solid lines are for the "relaxed" clusters while the 
dashed lines are for the "perturbed" clusters. Overlayed are the observational Suzaku point of George et al.(2009) for PKS0745-191 (diamonds), and of Bautz 
et al.(2009) for A1795 (triangles). 



scatter measured in relaxed systems. This is explained 
noticing that, even if clumps are efficiently removed by 
our filtering procedure, perturbed structures are gener- 
ally characterized by a larger amount of chaotic mo- 
tions and shocks at all radii (e.g. ValO for an analysis 
of shocks in the ENZO sample). An important point to 
notice here is that the agreement between the azimuthal 
scatter in two codes is generally better than the agree- 
ment between the simple radial profiles (Figj2]). This 
suggests that, despite of the small but non-negligible 
differences in the cluster profiles of the GADGET2 and 
ENZO samples, which can be due to a number of possi- 
ble reasons (i.e. slightly difference cosmological setup, 
cluster extractions procedures, numerical issues etc.), 
our predictions on the level of intrinsic azimuthal scat- 
ter in clusters are robust and independent on the partic- 
ular code adopted. 

It is now interesting to compare these findings with 
the (few) available observations of azimuthal scatter 
derived from long Suzaku exposures of PKS0745-191 
(George et al.2009, black diamonds in FigH]> and of 
A1795 (Bautz et al.2009, black triangles). The first 



observation was performed with Aa « ir/2 (thus 
N = 4), while for the second Aa « 7r (N = 2). 

In general, the reported azimuthal scatter of 
PKS0745-191 lies at the uppermost part of our aver- 
age distributions in FigfSJ while the azimuthal scatter 
of A 1795 lies in the lowermost part. Even if the radial 
sampling available to observations is yet quite scarce, it 
seems that the trends with radius of the available data 
can be reproduced by cosmological simulations, with 
the exception of the peculiar decreasing trend with ra- 
dius reported for the gas density of PKS0745-191. 

At least two reasons may be responsible for the 
larger azimuthal scatter of PKS0745-191 for r < -R200 
compared to our simulated results. First, PKS0745-191 
is a system affected by a un-usually large perturbation 
(caused for instance by a merger or by the accretion of 
a large scale filament), which is observed only in a few 
exceptional cases in our sample of simulated clusters. 
We notice, however, that PKS0745-191 is usually clas- 
sified as a relaxed cool-core cluster (George et al.2009 
and references therein); anyway this classification is 
just based on the cluster morphology around the core 
region, and thus might be not at variance with our defi- 
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Figure 6. Same as in Figure \5\ but for the un-filtered gas density and gas 
temperarure and for N = 8. 



nition of "perturbed" systems, for which the whole vol- 
ume of the cluster is analyzed. A second possibility, is 
that some un-resolved substructures are present in one 
of the 4 FOV targeted by Suzaku, a possiblity that the 
authors themselves considered in their work (George et 
al.2009). Indeed, if the azimuthal scatter is computed 
for the un-filtered gas density and Tsl (see Figj6], for 
all the simulations in our data-set, much larger scatters, 
compatible with the findings of George et al.(2009), are 
found. In the case of ENZO runs, moreover, it seems 
that on average also the "inverted" trend in gas density 
for 0.5 ^ r/i?200 ^ 1 is statistically reproduced by the 
simulated data, due to the effect of resolved filaments 
of gas matter impinging at the cluster periphery along 
some preferential sectors, which are well resolved with 
the mesh refinement strategy adopted in this particular 
data-set (see ValO for details). Similar results were also 
reported by Burns et al.(2010). 

This test compared the "intrinsic" azimuthal scat- 
ter measured in our filtered simulations (using a num- 
ber of sectors which produces a convergence in the es- 
timate of S c (r)) with the presently best available data 
from real Suzaku observations; however the data from 
George et al.(2009) were derived from N = 4 sectors, 
and the data from Bautz et al.(2009) were derived from 
N = 2. As we shown in Figf$l the use of different an- 
gular sampling is expected to affect the estimate of the 
real S c (r) for Aa > 7r/4. To investigate this effect, we 
also repeated the measures of the azimuthal scatter pro- 
files across the whole sample of simulations, assuming 
N = 2 (top panels of Fig|7]) and N = 4 (bottom pa nels 
of FigjT],. As expected from the tests shown in Sec |2.2| 
the azimuthal scatters are progressively reduced at all 
radii as the angular size of sectors is increased. With 
N = 2 our measurement of the average azimuthal scat- 
ter in clusters is now well in agrement with the data ob- 
tained for A1795 from Bautz et al.(2009) observations. 
In the case of = 4 still the scatter in gas density for 
PKS0745-191 from George et al.(2009), derived with 
the same number of sectors, can be explained by our 
data only assuming the contribution from un-resolved 
gas clumps. 



4 DISCUSSION AND CONCLUSIONS 

In this work we report upon the first statistical mea- 
sure of the azimuthal scatter in the thermal distribution 
of galaxy clusters (analyzing gas density, gas tempera- 
ture, gas entropy and X-ray emission) with a large set 
of cosmological simulations run with ENZO (Vazza et 
al.2010) and with GADGET2 (Roncarelli et al.2006). 
This joint sample of simulations offers a unique view 
to statistically constrain the intrinsic level of azimuthal 
scatter in the cluster observables, as produced by the 
simple gravitational processes which follow the hier- 
archical growth of cosmic structures. This can be com- 
pared to the (few) real data recently obtained with long 
Suzaku exposures of specific sectors in nearby clusters 
(e.g. PKS0745-191 by George et al.2009; A1795 by 
Bautz et al.2009). For this purpose, we apply a post- 
processing technique to our simulations designed to re- 
move the dense collapsed gas clumps which are usually 
masked out from real X-ray observations (Roncarelli 
et al.2006). According to our study, a statistically con- 
vergent modeling of the intrinsic azimuthal scatter in 
cluster profiles is achieved for sectors of angular size 
Aa « 7r/4 or smaller, almost regardless of the numer- 
ical scheme and of the dynamical state of clusters. The 
level of azimuthal scatter reported for the available ob- 
servations (George et al.2009; Bautz et al.2009) can be 
in principle well reconciled with our simulations when 
using a similar number of angular sectors (with angu- 
lar size Aa = vr/2 for PKS0745-191 and Aa = vr 
for A1795); however the peculiar trend of gas density 
for r < i?200 in PKS0745-191 may be well reconciled 
with simulations only in the case that some un-resolved 
gas substructure (e.g. a filament or a gas clump) is be 
present in one of the Suzaku FOVs. 

To summarize our results (for simplicity restricting 
to gas density and spectroscopic-like temperature), we 
report in Table 2 our estimates of the azimuthal scatters 
in our sample of clusters, after the "99 percent filter- 
ing" and assuming N = 8 and N = 4. 

In the future, a larger sample of real cluster ob- 
servations would be desirable in order to fully account 
for the variety of dynamical states present in clusters, 
which may affect the observed azimuthal scatter in a 
sizable ways. 

We also comment that, at the present stage, the 
effect of non-gravitational physical processes acting 
at the cluster peripheral regions (e.g. AGN-feedback, 
magnetic fields, Cosmic Rays, etc.) is not required to 
reproduce the presently available observations, given 
that present cosmological high-resolution and "simple" 
non-radiative runs are suitable to reproduce the trends 
reported by observations (see also the tests in the Ap- 
pendix). This findings are qualitative in line with pre- 
vious results discussed by Romeo et al.(2006); Roncar- 
elli et al.(2006); Burns et al.(2010). 

However, in the future a self-consistent modeling 
of the ICM metallicity and chemical enrichment pro- 
cesses in the cluster outskirts may be important to 
model the X-ray luminosity in a more realistic way 
(e.g. by adding chemical evolution, transport of met- 
als and line cooling). 
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Table 2. Reference values for the azimuthal scatter of gas density and spec- 
troscopic like temperature for the GADGET2 and ENZO samples, using the 
"99 percent filtering" at 3 different radius from the center, assuming N = 4 
and N = 8. The mean azimuthal scatter values are reported with their ±<r 
standard error. The last column shows the dynamical state of the cluster 
sub-samples considered (R=relaxed, P=perturbed). 
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Figure 7. Azimuthal scatter profiles for gas density and gas temperature for 
the same samples of Fig|5] for N=2 (top panels) and N=4 (bottom panels). 
The meaning of colors and linestyles is an in Fig|5] The points derived from 
A1795 (Bautz et al.2009) are reported in the top panels (as triangles), while 
the points derived from PKS0745-191 (George et al.2009) are reported in 
the bottom panels (as squares). 



ACKNOWLEDGEMENTS 

We thank C.Gheller and R.Brunino for their price- 
less support at HPC-CINECA (Italy), and G.Brunetti 
of fruitful discussions. We acknowledge the help of 
D.Fabjan in the treatment of GADGET simulations 
data. F.V. acknowledges the usage of computational 
resources under the CINECA-INAF 2008-2010 agree- 
ment and the 2009 Key Project "Turbulence, shocks 
and cosmic rays electrons in massive galaxy clusters 
at high resolution", and partial financial support from 
the grant FOR 1254 from Deutschen Forschungsge- 
meinschaft. S.E. and M.R. acknowledge the financial 
contribution from contracts ASI-INAF 1/023/05/0 and 
1/088/06/0. M. R. acknowledges the support of grant 
ANR-06-JCJC-0141. K.D. acknowledges the support 
by the DFG Priority Programme 1177 and additional 
support by the DFG Cluster of Excellence "Origin and 
Structure of the Universe". 



REFERENCES 

Bautz, M. W., et al. 2009, PASJ, 61, 1117 
Burns, J. O., Skillman, S. W, & O'Shea, B. W. 2010, 
ApJ, 721, 1105 

Cassano, R., Ettori, S., Giacintucci, S., Brunetti, G., 
Markevitch, M., Venturi, T, & Gitti, M. 2010, ApJL, 
721, L82 

Cavaliere , A., Lap i, A., & Fusco-Femiano, R. 2010, 
larXiv: 1010.44151 

Dolag K., Vazza F, Brunetti G., Tormen G., 2005, 
MNRAS, 364, 753 

Ettori, S., Fabian, A. C, & White, D. A. 1998, MN- 
RAS, 300, 837 

Ettori, S., & Balestra, I. 2009, A&A, 496, 3 43 
Ettori, S., & Molendi, S. 2010, arXiv: 1005.0382 
Frenk, C. S., et al. 1999, ApJ , 525, 554 
George, M. R., Fabian, A. C, Sanders, J. S., Young, 

A. J., & Russell, H. R. 2009, MNRAS, 395, 657 
Hoshino, A., et al. 2010, arXiv: 1001.5T33l 
Kawahara, H., Reese, E. D., Kitayama, T, Sasaki, S., 

& Suto, Y. 2008, ApJ, 687, 936 
Kawaharada, M., et al. 2010, ApJ, 714, 423 
Lapi, A., Fusco-Femiano, R., & Cavaliere, A. 2010, 

A&A, 516, A34 

Leccardi, A., & Molendi, S. 2008, A&A, 486, 359 
Liedahl, D. A., Osterheld, A. L., & Goldstein, W. H. 

1995, ApJL, 438, LI 15 
Loken, C, Norman, M. L., Nelson, E., Burns, J., 

Bryan, G. L., & Motl, P. 2002, ApJ, 579, 571 
Mitchell, N. L., McCarthy, I. G., Bower, R. G., The- 

uns, T, & Crain, R. A. 2009, MNRAS, 395, 180 
Neumann, D. M. 2005,A&A, 439, 465 
Norman M. L., Bryan G. L., Harkness R., Bordner 

J., Reynolds D., O'Shea B., Wagner R., 2007, arXiv, 

705, arX iv:0705.1 556 

Rasia, E., Mazzotta, P., Borgani, S., Moscardini, L., 
Dolag, K., Tormen, G., Diaferio, A., & Murante, G. 
2005, ApJL, 618, LI 

Romeo, A. D., Sommer-Larsen, J., Portinari, L., & 
Antonuccio-Delogu, V. 2006, MNRAS, 371, 548 



The scatter of clusters profiles. 9 



Roncarelli, M., Ettori, S., Dolag, K., Moscardini, L., 
Borgani, S., & Murante, G. 2006, MNRAS, 373, 
1339 

Reiprich, T. H., et al. 2009, A&A, 501, 899 

Springel, V., & Hernquist, L. 2003, MNRAS, 339, 289 

Springel V., 2005, MNRAS, 364, 1 105 

Springel V., 2010, MNRAS, 401, 791 

Valdarnini, R. 2010, arXiv : 1010.33781 

Vazza, F., Brunetti, G., Gheller, C, & Brunino, R. 
2010, New A, 15, 695 

Vazza, F. 2010, MNRAS, 1724 

Voit, G. M., Kay, S. T., & Bryan, G. L. 2005, MN- 
RAS, 364, 909 

Vikhlinin, A., Forman, W., & Jones, C. 1999, ApJ, 
525, 47 

Vikhlinin, A., Kravtsov, A., Forman, W., Jones, C, 
Markevitch, M., Murray, S. S., & Van Speybroeck, 
L. 2006, ApJ, 640, 691 

Wadsley, J. W., Veeravalli, G., & Couchman, H. M. R 
2008, MNRAS, 387, 427 



APPENDIX 

We briefly present here additional tests to assess the 
reliability of the results against a few numerical effects 
that might play a role in the computation of the average 
azimuthal scatter. 

First, we show the additional results for the pro- 
files and the azimuthal scatter for the the two most mas- 
sive clusters of the GADGET2 runs (gla and g8a), in 
re-simulations where cooling, star formation and feed- 
back from galactic winds where self-consistently fol- 
lowed during the run. The star formation is followed 
by adopting a sub-resolution multiphase model for the 
interstellar medium, including also the feedback from 
supernovae and galactic outflows (Springel & Hern- 
quist 2003). In these runs, the efficiency of supernovae 
to power galactic winds has been set to 50 per cent, 
which turns into a wind speed of ~ 340km/ s. These 
re-simulations have been extracted from the same "cfs" 
dataset of Roncarelli et al.(2006). In the top panel 
of FigfTJ we show the average gas density (top lines) 
and gas entropy (bottom lines) profiles for gla and 
g8a at z- 0, after the "99 percent" filtering procedure 
(Sec |2.2| ), for the reference physical run adopted in 
the main article (solid lines) and for the csf model 
(dashed). Except for the sizable increase of the core 
gas density and the decrease of the core gas entropy, as 
an obvious effect of radiative cooling, the agreement 
for all the radii outside r > 0.2i?2oo is remarkable. 
This is in line with the early results already reported by 
Roncarelli et al.(2006) and Romeo et al.(2006). 

In the bottom panel of FigfJJ we show the corre- 
sponding radial trend of the azimuthal scatter as de- 
fined in EqfJJ and keeping N = 8 sectors (the bold 
lines refer to gas density and the thin lines refer to 
gas entropy). For both clusters in this case the dif- 
ferences in the azimuthal scatter are sizable when the 
two physical runs are compared (~ 10 per cent at all 




Figure 1. Left panel: radial profiles for gas density (top lines) and gas en- 
tropy (bottom lines) for clusters gla and g8a of the GADGET2 sample. The 
solid lines refer to non-radiative runs, while the dashed lines refer to runs 
with radiative cooling, star formation and feedback from galactic winds. 
Right panel: profile of azimuthal scatter for N = 8 sectors. The thick lines 
refer to gas density, the thin lines refer to gas entropy; the meaning of colors 
and line-styles is as in the Left panel. 



radii); however all profiles scatter around the value of 
S(r) ~ 0.10 - 0.15 at all radii from r > 0.5R 20 o, 
with no evident trend with radius. This may be well ex- 
plained by the small scale effect of cooling sub-clumps 
at all radii in the csf run: the "99 per cent filtering" pro- 
cedure outlined in the main paper efficiently removes 
the pointlike contribution from these clumps, however 
cooling make them more compact and enhances their 
survival within the cluster atmosphere and their effi- 
ciency in injecting shocks and small scale disturbances 
in the ICM (e.g. Valdarnini 2010), which can slightly 
change the local scatter behaviour, even if not the over- 
all trend of azimuthal scatter in the clusters. These re- 
sults make us confident that the adoption of more com- 
plex physical modeling compared to the non-radiative 
runs investigated in the main paper cannot change our 
conclusion in any sizable way, for all the cluster vol- 
ume outside of cores. 

In a second test, we checked the stability of our 
results against the adoption of a larger line of sight 
across the projected cylinder, increasing its size from 
4 • -R200 ( as use d in the main paper) to 6 • R200 ■ This was 
motivated by the fact that the overlap of not virialized 
regions and large scale environment along the line of 
sight towards the cluster may vary the reported values 
of the azimuthal scatter. We thus re-computed the az- 
imuthal scatter of gas density, gas temperature and gas 
entropy for cluster gla of the GADGET2 sample, by 
considering a line of sight of 6 • -R200 (FigJU). The test 
clearly shows that the difference in the azimuthal scat- 
ter of all quantities remains essentially the same within 
a few percent, with all trends un-altered. We thus con- 
clude that our results are well constrained for a total 
line of sight of ^ 4 • i?200- 
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Figure 2. Radial profile of azimuthal scatter of gas density, gas temperature 
and gas entropy for cluster gla (JV = 8 sectors). The solid lines show the 
results for a line of sight of 4 ■ R200 . the dashed lines show the results for a 
line of sight of 6 ■ R200 ■ 



